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NOTATION 

Coefficient  in  general  auxiliary  equation,. 

Equation  [84] 

Constant  in  Equation  [19] 

Coefficient  in  general  auxiliary  equation, 

Equation  [84] 

Constant  in  Equation  [T9] 

Constant 

Drag  coefficient,  D/|-pU2S 
Drag 

Function  of  H  in  Table  3 

Base  of  natural  logarithms 

Shape  parameter  of  velocity  profiles,  6*/9 

Shearing-stress  distribution  integral, 

Subscript  referring,  to  Initial 
conditions 

Expcnent  of  power  law  for  velocity  profiles, 
Equation  [8] 

Subscript  for  flat-plate  quantities 
Pressure  in  fluid 
Exponent  in  Equation  [62] 

Reynolds  number  based  on  x,  Ux/v 
Reynolds  number  based  on  6,  U 6/v 
Transverse  radius  to  surface  of  body  of  revolution 
Surface  area 

Constant  in  Equation  [26] 

Subscript  indicating  quantities  at  the  trailing 
edge  of  an  airfoil 

Velocity  at  outer  edge  of  boundary  layer 

Undisturbed  velocity  of  fluid 

Velocity  in  boundary  layer  parallel  to  surface 

Friction  velocity,  frJ~P 

Velocity  in  boundary  layer  normal  tc  surface 

Width 

Coordinate  parallel  to  surface 
Coordinate  normal  to  surface 
Pressure-gradient  parameter,  — 

rw .  , 

Shearing-stress  parameter,  — 

T*  6i 

Burl  pressure-gradient  parameter,  -g  r)/4 


1 

I 


ill 


y  Velocity  parameter, 

y*  * 

6  Thickness  of  boundary  layer 

6 *  Displacement  thickness  of  boundary  layer,  J  (l-jj)dy 

(  Burl  skin-friction  parameter 

t)  Gruschwltz  velocity  profile  shape  parameter,  1  -  y* 

9  Momentum  thickness  of  boundary  layer,  J4jj(l-jj-)dy 

v  Kinematic  viscosity  of  fluid 

p  Density  of  fluid 

r  Shearing  stress 

rw  Shearing  stress  at  wall 

^  Parameter  for  integral  of  shearing-stresB 

distribution,/^4  ^  d  (f.) 
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ABSTRACT 

A  method  is  described  for  the  calculation  of  the  development  of  incompress¬ 
ible  turbulent  boundary  layers  in  pressure  gradients,  so  that  drag  and  separation  of 
flow  may  be  predicted  more  accurately.  The  principal  contribution  is  the  formulation 
of  new  relations  for  the  variation  with  pressure  gradient  of  1,  the  shearing  stress  at 
the  wall  (local  skin  friction)  and  of  2,  the  Integral  of  the  shearing  stresses  across 
the  boundary  layer  for  use  in  a  moment-of -momentum  equation  recently  derived  by 
Tetervln  and  Lin. 

With  the  assumption  of  the  one-parameter  characterization  of  velocity  profiles 
by  a  shape  parameter  H,  the  moment-of- momentum  equation  and  the  von  Karman  mo¬ 
mentum  equation  constitute  a  working  method  involving  the  solution  of  a  pair  of  simul¬ 
taneous  first-order  differential  equations. 

As  an  introduction  to  the  subject  of  turbulent  boundary  layers  In  pressure 
gradients,  a  review  of  previous  methods  Is  Included  in  this  report. 

INTRODUCTION 

An  accurate  calculation  of  the  growth  of  turbulent  boundary  layers 
in  a  pressure  gradient  is  required  for  many  problems  in  hydrodynamics,  aero¬ 
dynamics  and  hydraulics  where  high  Reynolds  numbers  are  involved.  Some  of 
these  problems  are  the  calculation  of  the  drag  of  airfoils  and  of  bodies  of 
revolution,  the  prediction  of  separation  of  flow  on  airfoils,  the  prediction 
of  flow  conditions  in  the  test  sections  and  diffusers  of  water  tunnels  and 
wind  tunnels,  and  the  prediction  of  the  flow  (frictional  wake)  around  the 
tails  of  bodies  of  revolution  for  use  in  designing  propellers  and  control 
surfaces. 

Unlike  laminar  boundary  layers,  turbulent  boundary  layers  are  not 
very  amenable  to  mathematical  analysis  owing  to  the  present  incomplete  stat*» 
of  knowledge  of  the  mechanics  of  turbulent  flow.  Nevertheless,  considerable 
success  in  analyzing  turbulent  boundary  layers  on  flat  plates  with  zero  pres- 
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sure  gradient  has  been  achieved  by  such  investigators  as  Prandtl  and  von 
Karman.  No  comparable  progress  has  been  attained  to  date,  however,  in  the 
theory  of  turbulent  boundary  layers  in  a  pressure  gradient,  especially  in  an 
adverse  (positive)  pressure  gradient. 

Briefly,  the  current  objective  for  analyzing  turbulent  boundary 
layers  in  a  pressure  gradient  is  to  develop  an  equation  for  the  variation  of 
the  shape  of  the  velocity  profile  within  the  boundary  layer.  This  formulation 
is  required  as  an  auxiliary  equation  for  complementing  the  well-known 
von  Karman  momentum  equation.  Auxiliary  equations  have  been  obtained  empir¬ 
ically  by  a  number  of  Investigators.  Because  of  significant  differences  ex¬ 
isting  among  these  empirical  formulations,  Tetervin  and  Lin  attempted  a  the¬ 
oretical  approach  to  the  problem.  One  of  the  equations  obtained  by  them  is'  a 
moment-of -momentum  equation  derived  from  the  basic  Prandtl  boundary-layer 
equations  and  from  a  one-parameter  characterization  for  the  shape  of  the  ve¬ 
locity  profiles.  Although  the  moment-of -momentum  equation  of  Tetervin  and  Lin 
has  the  required  form  for  an  auxiliary  equation,  it  lacks  essential  relations 
needed  for  numerical  computations— namely,  the  variation  with  pressure  gra¬ 
dient  of  the  shearing  stress  at  the  surface  (local  skin  friction)  and  of  the 
integral  of  the  distribution  of  the  shearing  stresses  across  the  boundary 
layer.  The  principal  aim  of  this  paper  is  to  supply  the  necessary  shearing- 
stress  relations  to  this  moment -of -momentum  equation  in  order  to  develop  a 
suitable  auxiliary  equation  with  a  theoretical  basis,  to  be  used  In  conjunc¬ 
tion  with  the  von  Karman  momentum  equation. 

The  relation  for  the  local  skin  friction  in  a  pressure  gradient  Is 
derived  on  the  basis  of  the  recent  work  of  Ludwieg  and  Tlllmann  who  demon¬ 
strated  the  validity  of  applying  the  so-called  "law  of  the  wall"  to  the  inner 
flow  in  the  boundary  layer.  The  integral  of  the  transverse  distribution  of 
the  shearing  stresses  in  a  pressure  gradient  is  found  empirically  in  this 
paper  to  be  a  function  of  similar  flat-plate  data  for  zero  pressure  gradient. 

The  moment -of -momentum  equation,  after  being  modified  by  the  inclu¬ 
sion  of  the  shearing-stress  relations,  is  demonstrated  to  agree  with  the  aver¬ 
age  of  other  existing  auxiliary  equations  whose  formulation  is  based  on  purely 
empirical  grounds. 

FUNDAMENTAL  RELATIONS  FOR  TURBULENT  BOUNDARY  LAYERS 

Some  of  the  fundamentals  of  turbulent  boundary  layers  are  briefly 
reviewed  here  as  an  introduction  to  the  subject.  The  comprehensive  summaries 
by  Prandtl,  Ooldstein,  and  Schllchting,  References  1,  2,  and  5  respectively, 
should  be  referred  to  for  fuller  details.* 
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In  general,  steady  fluid  flow  (quasi-steady  for  turbulent  flow), 
manifests  Itself  In  two  distinct  types  of  motion,  laminar  at  low  Reynolds  num¬ 
bers  and  turbulent  at  high  Reynolds  numbers.  Whereas  laminar  flow  proceeds 
in  a  regular  pattern  of  streamlines,  turbulent  flow  advances  In  a  haphazard 
combination  of  mixing  motions.  Although  the  variations  in  velocity  and  pres¬ 
sure  of  turbulent  flow  follow  a  random  course,  they  may  be  considered  to  fluc¬ 
tuate  about  some  mean  value  as  shown  in  Figure  1 . 

The  analytic  treatment  of  turbulent  flow  consists  in  separating  the 
fluid  motion  into  a  mean  flow  and  a  superposed  fluctuation  flow.  Then  the 
mean  and  fluctuation  velocities  and 
pressures  are  substituted  into  the 
Navler-Stokes  equations  for  viscous 
flow  and  into  the  equation  of  conti¬ 
nuity.  After  appropriate  time  aver¬ 
ages  are  formed,  the  resulting  equa¬ 
tions  in  terms  of  mean  quantities 
resemble  those  for  laminar  flow,  with 
the  exception  of  additional  terms  in¬ 
volving  averages  of  various  products 
of  the  fluctuation  quantities.  These 
fluctuation  terms  represent  the  mixing 
motion  of  turbulence  and  act  as  appar¬ 
ent  stresses  (Reynolds  stresses}  with¬ 
in  the  fluid.  Since  an  analytical 
description  of  the  Reynolds  stresses 
has  not  been  formulated,  a  direct  solu¬ 
tion  of  the  differential  equations  for 
turbulent  flow  is  virtually  impossible. 

The  theory  of  the  flow  of  viscous  fluids  serves  as  a  useful  guide,  however, 
even  though  empirical  or  semi-empirical  methods  may  be  followed  in  many  prac¬ 
tical  cases. 

At  relatively  high  Reynolds  numbers,  the  viscous  effects  of  fluids 
like  air  or  water  are  confined  mostly  to  a  narrow  region  or  bdnd  next  to  the 
solid  surface.  The  laminar  flow  through  this  region,  or  boundary  layer  as  it 
is  better  known,  acts  in  the  same  manner  as  fully  developed  viscous  flow  in 
pipes,  inasmuch  as  it  becomes  turbulent  after  a  critical  value  of  some  Reyn¬ 
olds  number  is  exceeded.  The  equations  of  motion  for  the  turbulent  boundary 
layer  are  found  analytically  from  the  Navier-Stokes  equations  for  turbulent 
flow  in  a  way  analogous  to  that  for  the  laminar  boundary  layer— that  is,  by 


Tim* - fc- 

Figure  1  -  Velocity  Fluctuations  of 
Turbulent  Flow  at  a  Fixed  Position 
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rejecting  terms  involving  small  orders  of  magnitude.  For  steady  incompress¬ 
ible  flow  with  negligible  longitudinal  curvature,  the  turbulent-boundary-layer 
equations  of  motion  are  found  to  be 
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for  0  s  y  <  6  and  0  <  u  s  U.  Here  u  and  v  are  the  x-  and  y-components  of  the 
velocity  within  the  boundary  layer,  parallel  and  normal  to  the  surface  of  the 
body  respectively,  U  is  the  velocity  at  the  outer  edge  of  the  boundary  layer, 
d  is  the  thickness  of  the  boundary  layer,  p  is  the  density  of  the  fluid,  and 
r  is  the  total  shearing  stress  within  the  fluid.  It  18  to  be  noted  that  the 
velocities  and  pressures  Indicated  in  Equation  [l]  represent  time  averages  of 
the  turbulent  quantities  in  the  flow.  In  addition  to  two-dimensional  flow, 
Equation  [1]  is  also  applicable  to  axisymmetrlc  flow  when  the  boundary-layer 
thickness  d  Is  small  relative  to  the  transverse  radius  of  the  solid  boundary 

rw 

The  accompanying  equation  of  continuity  for  two-dimensional  flow  is 


du  A  dv  n 

57  +  57“  0 


[2a  j 


and  for  axisymmetrlc  flow  it  is 


for  d  «  rw. 

From  the  boundary-layer  equation  of  motion,  Equation  [1],  and  the 
equation  of  continuity,  Equation  [2],  the  von  Karman  momentum  equation  is  de¬ 
rived  by  integrating  each  term  of  [1]  with  respect  to  the  y-coordlnate.  The 
von  Karman  equation  represents  the  basic  working  equation  for  calculating  the 
growth  of  boundary  layers  in  the  presence  of  pressure  gradients.  For  purposes 
of  analysis  the  momentum  equation  is  best  expressed  in  terms  of  the  displace¬ 
ment  thickness  d#,  the  momentum  thickness  9  and  the  shape  parameter  H  which 


are  defined  as 


d*  “  J0  (1  "  ir)dy  (3] 
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In  terms  of  these  quantities  the  von  Ka'rman  momentum  equation  for  two- 
dimensional  flow  becomes 

and  for  axisymmetric  flow  it  becomes  for  6  «  rw 
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dx 
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(H  +  2) 


V  dU 
U  dx 


[6b  1 


where  rw  is  the  shearing  stress  at  the  wall. 

Physically  the  von  Karman  momentum  equation  is  Interpreted  as  repre¬ 
senting  the  rate  of  change  of  the  momentum  of  the  fluid  within  the  boundary 
layer  as  a  function  of  the  frictional  resistance  and  of  the  pressure  gradient. 
When  solving  a  particular  problem.  U,  the  velocity  at  the  outer  edge  of  the 
boundary  layer,  and  dU/dx,  representing  the  pressure  gradient,  are  given  by 
measurements  from  pressure  taps  on  the  body  or  by  potential -flow  calculations. 
Quantities  such  as  the  shape  parameter  H  and  the  local-skin-frictlon  coeffi¬ 
cient  Tw/oUa,  have  to  be  determined  by  other  means.  By  assuming  flat-plate 
values  for  H  and  rw//»U* ,  the  momentum  equation  can  be  integrated  to  give  an 
approximate  solution  which  in  many  caBSB  is  not  adequate  for  the  problem  in 
question.  The  more  exact  solution  depends  on  acquiring  auxiliary  relations 
for  the  variation  of  H  and  r^/pM2  in  a  pressure  gradient .  The  detailed  con¬ 
sideration  of  H  andrw/pU2  in  a  pressure  gradient  is  the  major  theme  of  this 
paper  and  will  be  discussed  fully  in  the  subsequent  sections. 

As  an  illustration  of  the  flow  idealized  by  boundary-layer  theory, 
Pigure  2  shows  the  flow  pattern  around  a  body  of  revolution.  Essentially, 
there  is  potential  flow  outside  the  boundary  layer  and  viscous  flow  within. 

The  viscous  boundary-layer  flow  is  laminar  from  the  nose  up  to  the  transition 
zone  where  it  develops  into  a  turbulent  flow.  Beyond  the  tail,  the  flow 
through  the  boundary  layer  merges  to  form  a  wake  trailing  behind  the  body. 

The  velocity  U  at  the  outer  edge  of  the  boundary  layer  is  that  of 
potential  streamline  and  varies  from  zero  at  the  stagnation  point  at  the  nose, 
through  a  range  of  values  along  the  body  to  the  free-stream  velocity  U«>  in 
the  wake  at  infinity.  Although  the  pressure  p  ac  the  outer  edge  of  the  bound¬ 
ary  layer  is  essentially  the  same  as  that  across  the  boundary  layer,  it  varies 
along  the  body  with  x  and  thuB  produces  a  pressure  gradient.  The  pressure  p 
and  the  velocity  U  of  the  potential  streamline  are  related  by  the  Bernoulli 
equation  as  follows 
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or 


p  +  jpU8 


constant  'i 
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A  negative  pressure  gradient  is  favorable  and  a  positive  pressure  gradient 
adverse  for  preventing  separation  of  flow. 

Within  the  boundary  layer  the  velocity  u  is  zero  at  the  surface, 
since  no  slippage  occurs,  and  increases  to  U  at  the  outer  edge  of  the  boundary 
layer.  This  variation  of  u  with  y,  or  the  velocity  profile,  is  very  Important 
in  the  study  of  boundary  layers.  In  the  case  of  turbulent  boundary  layers, 
the  following  power  law  variation  gives  an  excellent  approximation  for  most 
flow  conditions,  and  a  rough  approximation  for  strong  adverse  pressure 
gradients. 


where  the  value  of  the  exponent  n  depends  on  the  previous  development  of  flow 
in  the  boundary  layer. 


Figure  2  -  Typical  Boundary  Layer  Around  a  Body  of  Revolution 
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EXISTING  METHODS  FOR  CALCULATING  TURBULENT  BOUNDARY  LAYERS 
IN  A  PRESSURE  GRADIENT 

Interest  in  effects  of  pressure  gradient  on  turbulent  boundary  lay¬ 
ers  has  centered  largely  on  adverse  (positive)  pressure  gradients  causing  sep¬ 
aration  of  flow  such  as  occur  in  diverging  conduits  and  on  the  suction  Bide 
of  airfoils  at  large  angles  of  attack.  (A  review  of  the  literature  on  this 
subject  appears  in  Reference  4.)  Since  separation  of  flow  is  characterized 
by  a  reversal  in  the  direction  of  flow  at  the  surface,  the  velocity  profile 
of  the  boundary  layer  suffers  marked  changes  in  shape  as  the  flow  approaches 
the  point  of  separation,  as  seen  in  Figure  3<  Hence  the  study  of  flow  lead¬ 
ing  to  separation  must  be  directed  towards  the  determination  of  the  manner  in 
which  the  shape  of  the  velocity  profile  is  affected  by  pressure  gradients,  as 
well  as  by  other  factors,  such  as  Reynolds  number. 

Since  a  single  parameter  for  the  shape  of  the  velocity  profile 
greatly  simplifies  the  mathematical  analysis,  considerable  effort  towards  this 
goal  has  been  expended  by  various  investigators.  Oruschwitz9  plotted  u/U 
against  y/o  in  his  experimental  study  of  velocity  profiles  and,  as  Bhown  in 
Figure  3,  obtained  a  family  of  velocity  profiles  in  various  pressure  gradients. 
He  concluded  that  the  value  of  u/U  at  Borne  standard  value  of  y/e  could  be  used 
to  characterize  the  shape  of  the  velocity  profile  for  any  pressure  gradient. 
Following  this  concept  Oruschwitz  defined  a  form  parameter  17  such  that 


V  •  1 


[91 


where  y  ■ 


Since  the  shape  parameter  H,  which  can  be  written 

-re-WD  1,01 


la  also  a  function  of  the  shape  of  the  velocity  profile,  it  cannot  be  inde¬ 
pendent  of  17  if  the  single-parameter  relationship  holds. 

This  is  confirmed  by  Oruschwitz  who  plotted  H  against  >7  for  a  series 
of  test  data.  As  seen  in  Figure  4,  a  fairly  smooth  variation  between  17  and  H 
results.  It  is  to  be  noted  that  Kehl,'  from  tests  covering  a  wider  range  of 
Reynolds  numbers  than  Oruschwitz ,  found  that  rj  is  not  a  unique  function  of  H 
but  that  it  varies  slightly  with  Reynolds  number. 


8 


Figure  3  -  Velocity  Profiles  of  a  Turbulent  Boundary  Layer 
In  a  Pressure  Gradient  (Gruschwitz,  Reference  5) 
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Figure  4  -  Shape  Parameter  H  as  Function  of  Qruschwitz  Form  Parameter  17 

(Schlichting,  Reference  3) 

Instead  of  rj  ,  von  Doenhoff  and  Tetervin7  preferred  the  use  of  the 
shape  parameter  H  Itself  as  the  single  parameter  to  characterize  the  shape  of  p 
the  velocity  profiles.  Points  from  an  extensive  collection  of  test  data  were 
plotted  together  in  the  form 

f-  r(f «)  i”) 


As  seen  In  Figure  5.  plots  of  u/U  against  H  yield  single  curves  for  each  value 
of  y/tf. 


Figure  5  -  Comparison  of  Power-Law  Velocity  Profiles  with 
Experimental  Points  of  von  Doenhoff  and  Tetervin, 
Reference  7 


For  purposes  of  mathematical  analysis  it  is  found  very  useful  to . 
consider  a  general-power-law  distribution  for  the  velocity  profiles  of  the 
form 

ir  •  (i)'  i»i 

tJslng  this  power  law  Pretsch8  derived  the  following  relations  between  the 
boundary- layer  parameters 


i*  n  H-l 

i  "  n+T  "  B+T 


[12] 


[14] 


H  -  2n  +  1 
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and  finally 

v  ‘  1  •  [hTSTii]  fl5j 


The  curve  for  the  relation  just  derived  for  r?(H)  Is  drawn  on  Figure  4  and 
shows  close  agreement  with  the  test  data  of  Qruschwitz.  Likewise  the  curve 
for  the  power  law  for  the  velocity  profiles  may  be  expressed  in  the  form 


[16] 


Curves  drawn  from  this  equation  in  Figure  5  show  close  agreement  with  the  data 
of  von  Doenhoff  and  Tetervin.  From  the  experimental  evidence  just  considered, 
it  is  seen  that  the  general  power  law  is  an  excellent  approximation  for  the 
velocity  profiles  of  turbulent  boundary  layers. 

After  the  particular  parameter  for  the  shape  of  the  velocity  profile 
has  been  selected,  its  variation  in  a  pressure  gradient  remains  to  be  deter¬ 
mined.  It  is  now  recognized  that  the  rate  of  variation  of  this  shape  param¬ 
eter  and  not  the  parameter  Itself  is  dependent  upon  the  local-pressure- 
gradient  parameter  and  upon  other  local  parameters  of  the  boundary  layer. 
However,  an  early  study  by  Burl9  in  1931  on  accelerated  and  retarded  flows  in 
closed  conduits  considered  the  shape  of  the  velocity  profile  to  depend  direct¬ 
ly  on  the  pressure  gradient. 

The  parameter  used  by  Burl  fbr  the  shape  of  the  velocity  profile  is 


r  m  —  — 
r  U  dx 


(17) 


where  R  «=  —  and  v  is  the  kinematic  viscosity  of  the  fluid.  Actually  r  Is  a 

•  V 

local  parameter  for  the  pressure  gradient  and  hence  can  only  roughly  represent 
the  shape  of  the  velocity  profile  whose  development  depends  on  the  previous 
history  of  the  flow.  The  parameter  r  as  used  by  Burl  in  his  analysis  has 
merit  in  providing  a  simplified  method  for  calculating  the  momentum  thickness 
0  approximately.  When  r  is  substituted  into  the  von  Karman  momentum  equation, 
Equation  [6a],  the  following  expression  is  obtained 


where  {  •  — \  r)/4  • 
PU‘ 

expression  to 


8x«r*i  -£[<-r("*})]  msi 

From  experimental  evidence,  Burl  linearized  this  ex-* 


(!»> 


•  ...  ,,tv  ...  .  , 

1 1 

or 

<L  (  r  1/4  a  \  _  a  dU/R 
dxlR«  9)  ~  a  Tf  dPR«  9) 

where  a  and  b  are  empirical  constants.  Equation  [19]  is  integrated  as  a 
linear  first-order  equation  in  terms  of  R ^*6  to  provide  a  direct  solution  for 
the  momentum  thickness  6.  Prom  test  data  Burl  considered  separation  to  occur 
at  r  -  -0.06.  The  Burl  method  of  analysis  was  continued  by  Howarth10  who  re¬ 
fined  the  method  by  solving  Equation  [  1 8 ]  directly,  using  experimental  data 
for  t  and  H. 

At  the  same  time  that  Buri  published  his  work,  Gruschwltz5  presented 

a  more  satisfactory  method  of  analysis.  Introducing  the  shape  parameter  r\ 

and  realizing  that  its  rate  of  change  in  a  pressure  gradient  is  more  important 

than  its  local  value,  Gruschwltz  plotted  the  nondlmensional  differential 

quantity  4—  ( U2 rj )  against  tj  and  obtained  a  linear  variation  for  his  test 
U  r 

data  as  shown  in  Figure  6.  Although  Gruschwltz  expected  some  effect  due  to 
Reynolds  number,  the  narrow  range  of  his  data,  1  x  103  <  R#  <  4  x  103 ,  pre¬ 
cluded  any  such  indication  and  consequently  he  included  no  Reynolds  number 


0.3  0.4  0.3  0.6  0.7  0.8  0.9  1.0 

V 

Figure  6  -  Gruschwltz  Form  Parameter  ij  as  Function  of 
Pressure  Gradient  (Gruschwltz,  Reference  5) 
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effect  in  his  analysis.  The  empirical  equation  fitting  his  test  points  is 
given  as 

_±  d(Uliz)  .  0 . 00894 -  0.00461  [20] 

u*  ax 


for  77  <  0.8.  When  —•  is  expanded  to  2??(-^  it  is  seen  that 

the  rate  of  variation  of  77,  drj/dx,  and  the  nondimenslonal  pressure  gradient 
"u  dx  are  ^nc^u<^e<^  *n  Equation  [20].  The  variation  of  77  is  calculated 
from  the  auxiliary  equation,  [20],  when  used  simultaneously  with  the 
von  Karman  momentum  equation.  Qruschwitz  considered  separation  to  occur  at 
1 1  ^  0.8  which  is  equivalent  to  H  ^  1 .85. 

From  test  data  covering  a  wider  range  of  Reynolds  numbers, 

1  x  10s  <  R(  <  3.5  *  104,  Kehl*  Incorporated  the  effect  of  Reynolds  number  in 
the  Qruschwitz  equation  and  obtained 


±  d(U2w) 
U2  “ 


0.00894»7 


0.0164  0.85 

■  i°g10V  TV300T 


[21] 


Using  H  rather  than  77  as  the.  single  parameter  for  the  shape  of  the 
velocity  profiles,  von  Doenhoff  and  Tetervln7  made  an  empirical  study  of  the 
variation  of  H  in  pressure  gradients  in  order  to  obtain  an  auxiliary  equation. 
From  a  collection  of  NACA  and  other  test  data  representing  a  range  of  Reynolds 
numbers, 1  x  10s  <  R#  <  7  x  104  roughly,  they  considered  the  rate  of  change  of 
H  to  vary  as  follows 


•§-  [-*•«»»  -  '-286)  -ttfSp*0"-1-””  WJ 

po* 

where  r  /pU2  is  taken  from  the  flat  plate  formula  of  Squire  and  Young11 
"0 

jjja  *  [5.890  log',  (4.075  R,)]  [23] 

Von  Doenhoff  and  Tetervln  considered  separation  to  occur  in  the  range 
1.8  <  H  <  2.6.  It  is  to  be  noted  that  the  pressure-gradient  parameter  is 

and  the  parameter  representing  the  Reynolds-number  effect  is  rw  /plf2,  the 
local -skin-friction  coefficient.  Although  these  parameters  are  satisfactory 
in  themselves  for  describing  the  flow,  their  combination  into  a  single  term  in 
the  empirical  equation  seems  objectionable  for  analysis. 


15 


Using  substantially  the  same  data  as  von  Doenhoff  and  Tetervin  and 
a  somewhat  different  analysis,  Garner12  obtained  a  more  amenable  auxiliary 
equation  for  the  rate  of  variation  of  H 

-.dH  r  0.0T35(H-1.4)  •dUl.5(H-i.4) 

®dx“  [ - -  H  d3Eje 

The  following  flat-plate  formula  of  Falkner13  for  the  skin-friction- 
eoeffleient  la  used  in  ita  determination 


0.006534 


[25] 


A  noteworthy  feature  of  Garner's  equation  is  the  separation  of  the  pressure- 
gradient  parameter  -g-  ^  and  the  Reynolds  number  R#  into  different  terms. 

Theoretical  derivations  of  an  auxiliary  equation  have  proved  very 
difficult.  A  strenuous  but  unsuccessful  effort  was  made  by  Coleman14  in  1 9^+7 - 
A  very  promising  avenue  of  approach  was  introduced,  however,  more  recently  by 
Tetervin  and  Lin15  who  Investigated  various  integral  forms  of  the  Prandtl 
boundary-layer  equation,  [1).  By  multiplying  each  term  of  this  equation  by 
an  arbitrary  function  and  then  integrating,  it  is  possible  to  derive  various 
families  of  equations.  Using  as  the  multiplying  function  uSl  ys* ,  where  sx 
and  s2  are  arbitrary  constants,  the  general  integral  form  of  the  boundary- 
layer  equation  is  written 


fuSx  u8‘y8*dy  +  J\  &  uV’dy  «  -I  jf  &  u8>  y8*dy  +  usiys’dy 


[26] 


As  before,  this  equation  is  accompanied  by  the  equation  of  continuity,  [2]. 
When  both  exponents  sx  and  s2  are  set  equal  to  zero,  the  von  Karman  momentum 
equation  is  specified.  When  sx  ■  1  and  s2  *  0,  an  energy  equation  results; 
when  sx  =  0  and  s2  «  1 ,  a  moment- of-momentum  equation  is  indicated. 

The  moment-of-momentum  equation  appears  to  be  the  most  promising  of 
•these  derived  integrated  equations  for  obtaining  a  theoretical  auxiliary  equa 
tlon  for  the  variation  of  H.  Assuming  the  general-power-law  distribution  for 
the  velocity  profile,  [8],  Tetervin  and  Lin  obtained  the  moment -of -momentum 
equation  in  the  form 


-  dH  H(H+1)(H*-1)  9  dU  .  u/uS 
9  dT“ - ■  TJHx  + 


uV 

pu2 


(H  +  1)(Ha  -  l)fl-^-d(i)  [27] 

J0pU2 
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This  equation  is  applicable  to  two-dimensional  flows  and  to  axl- 
symmetric  flows  (d  «  rw) .  For  purposes  of  calculation  it  is  seen  that  rela¬ 
tions  for  the  variation  of  the  shearing  stresses  in  a  pressure  gradient  are 
still  required.  The  main  objective  of  this  paper  is  the  conversion  of  the 
above  equation  into  a  usable  auxiliary  equation  by  supplying  relations  for  the 
variation  in  a  pressure  gradient  of  the  local-skinrfrictlon  coefficient  r^/pU* 
and  the  integral  of  the  shearing -stress  distribution  f-^d^J. 

For  comparison  the  various  auxiliary  equations  for  the  variation  of 
the  shape  parameter  in  a  pressure  gradient  are  listed  in  Table  1 .  It  is  seen 
that  Burl's  method  does  not  provide  for  an  auxiliary  equation  since  no  allow¬ 
ance  is  made  for  the  effect  of  the  previous  development  of  the  flow.  On  the 
other  hsnd,  the  Gruschwitz  method  is  correct  in  utilizing  the  rate  of  change 
of  a  shape  parameter  in  its  analytic  formulation.  It  is  seen  that  Kehl  re¬ 
moved  one  of  the  main  objections  to  the  Gruschwitz  method  by  incorporating  the 


TABLE  1 

Summary  of  Methods  for  Calculating  the  Turbulent  Boundary  Layer 

in  a  Pressure  Gradient 


Investigator 

Criterion 

for 

Separation 

Auxiliary  Equation 

(for  uaa  Kith  von  Kdradn  Honantun  Equation) 

Raaarka 

Burl.  1931 
Reference  9 

r-  -0.06 

Iona  explicitly;  von  Karnan  aonantun  aquation  aodlfled  to 

E  (R?  *  )  ■  •  -  T  [19) 

r  la  not  proper  param¬ 
eter  for  shape  of  ve¬ 
locity  profile;  no 
allowance  in  aethod 
for  affect  of  previous 
development  of  flow  in 
boundary  layer 

Oruechwlta, 

1931 

Rifvaxi  J 

a  >0.8 

~L  ligli  .  0.00894a  -  0.00461  (20) 

lr 

No  Reynolda-nuaber  af¬ 
fect  Included;  baaed 
on  teat  data, 

1  x  10*<  R.  <  4  x  10* 

Haiti.  1943 

Rafaranea  6 

0.00894,.  (21) 

mu 

von  Doanhoff 
and  Tatarvln, 
1943, 

Rafaranea  7 

l.8<H<2.6 

[-2.039(8  :  1.286)  (^^Ja^-Mia)  [22j 

^  '[3.890  log*’ (4.079  «,)]  f  J 

kplrleal;  based  on 
taat  data, 

1  x  10*<R#<7  *  10* 

Oarnar,  19**** 
Rafaranea  12 

IBS 

Baaed  aottly  on  sane 
data  as  von  Doanhoff 
and  Tatarvln 

Tatarvln  and 
Un,  1990. 
Rafaranea  19 

Theoretically  'derived 
but  lacking  in  shear¬ 
ing-stress  lnforaatlon 
required  for  nunarleal 
solution 

(ihii  n  »»«**■>•:  f.v*r~-*»**  o>-». 


vwm*maa#\imm  ^we^tw 
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effects  of  Reynolds  number.  Mention  is  made  of  an  experimental  study  of  the 
Gruschwltz  method  by  Peters16  who  found  that  it  gave  fair  results  for  rj  but 
that  it  seemed  useless  for  predicting  separation.  The  form  of  the  Garner 
equation  is  superior  to  that  of  the  von  Doenhoff  and  Tetervin  equations  be¬ 
cause  it  separates  the  pressure-gradient  parameter  and  the  Reynolds -number 
effect.  A  critical  comparison  of  the  two  equations  just  mentioned,  made  by 
Fage  and  Raymer,17  claimed  serious  numerical  discrepancies  between  them  even 
though  they  had  both  been  derived  for  the  most  part  from  the  same  test  data . 
The  moment -of -momentum  equation  developed  in  the  present  report  will  be  found 
similar  in  form  to  that  of  the  Garner  equation.  A  graphical  comparison  of  all 
the  methods  with  the  exception  of  the  Burl  method  is  shown  at  the  end  of  this 
report . 


LOCAL  SKIN  FRICTION  AS  A  FUNCTION  OF  PRESSURE  GRADIENT 

The  manner  of  the  variation  of  the  local  skin  friction  or  shearing 
stress  at  the  wall  rw  with  pressure  gradient  as  well  as  with  Reynolds  number 
is  required  in  both  the  von  Karmen  momentum  equation  and  the  moment-of- 
momentum  equation. ^  Until  recently  the  experimental  study  of  wall  shearing 
stresses  was  hindered  by  the  poor  precision  of  the  existing  experimental  tech¬ 
niques,  based  on  measuring  the  rate  of  change  of  the  momentum  and  the  pressure 
of  the  flow  through  the  boundary  layer.  In  19^9  Ludwleg  and  Tlllmann18  intro¬ 
duced  a  new  method  for  accurately  measuring  the  wall  shearing  stress  in  a  pres¬ 
sure  gradient.  The  improved  technique,  described  in  Reference  19.  consists  in 
measuring  the  rate  of  heat  transferred  to  the  fluid  from  a  calibrated  instru¬ 
ment  Imbedded  in  the  wall.  From  their  measurements  Ludwieg  and  Tlllmann 
proved  that  the  wall  shearing  stress  diminishes  in  an  adverse  pressure  gradi¬ 
ent,  even  close  to  the  separation  point.  In  addition,  they  demonstrated  that 
the  so-called  "law  of  the  wall"  can  be  used  for  a  quantitative  analysis  of  the 
variation  of  the  wall  shearing  stress  in  a  pressure  gradient. 

The  "law  of  the  wall"20  states  that  within  the  boundary  layer  in  a 
pressure  gradient  there  is  a  region  of  flow  next  to  the  wall  which  has  veloc¬ 
ity  and  skin-friction  characteristics  similar  to  that  of  a  flat  plate  in  a 
zero  pressure  gradient.  As  seen  in  Figure  7  this  region  of  flow  or  inner  flow 
has  subregions  of  a  laminar  sublayer  and  a  transition  zone.  The  Prandtl 
friction-velocity  relationship  for  a  flat  plate2  can  then  be  applied  to  the 
whole  inner  flow  or 


[28] 


i6 


where  uT  «|/^~ls  the  friction  velocity.  Tne  functional  relationship  satisfy¬ 
ing  the.  turbulent  portion  of  the  inner  flow  Is 


Ci  +  c8 


[29) 


For  simplicity  of  analysis,  the  above  logarithmic  formulation  can  be  closely 
approximated  by. the  general  power  law 


uTyx"o 


■-m 


[30] 


By  virtue  of  the  law  of  the  wall,  C  and  n#  are  independent  of  pressure  gradi¬ 
ent  and  depend  only  on  the  Reynolds  number.  The  applicability  of  the  power- 
law  formulation  is  demonstrated  in  Figure  8  where  a  logarithmic  plot  of  u/U 
against  y/o  is  shown  for  various  pressure  gradients.  It  is  seen  that  the 
curves  are  straight  and  parallel  up  to  about  y  ■  B  .  This  situation  can  be  ex¬ 
plained  analytically  from  the  power  law  as  follows.  Substituting  y  ■  e  and 
u  ■  u#  into  the  power  law  gives 


Ha 

uT 


Dividing  Equation  [30]  by  [31]  produces 


or 


[31] 


[32] 

[33] 


Since  n#  does  not  depend  on  pressure  gradient,  the  curves  of  log  u/U  should  be 
straight  and  parallel  up  to  y  -  9. 

Ludwieg  and  Tlllmann  demonstrated  experimentally  that  the  law  of  the 
wall  is  also  valid  for  analyzing  the  variation  of  the  wall  shearing-  stress  in 
a  pressure  gradient.  From  test  data  in  the  range  of  Reynolds  number 
1  x  109  4  R#  <  4  x  104,  they  obtained  the  following  expression  for  the  varia¬ 
tion  of  the  local-skin-friction  coefficient  in  a  pressure  gradient 


Zl  s  0-0290  y1.705 
pU*  R°’*8 


[34] 


-0.398H 


where 


y  =  2.533  x  10 
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A  general  equation  for  the  variation  of  the  skin  iriction  in  a  pres¬ 
sure  gradient  will  now  be  derived  from  the  law  of  the  wall  by  using  the  power  - 
law  formula  for  the  velocity  profile  of -the  inner  flow, 


u_ 

“r 


"o 


[30] 


As  before,  at  y  ■  9  ,  u  ■  u#  and 


UT0x»O 


—  =  c(  — 1 — ) 

uT  \  v  I 


After  substituting  ur 


|/^"  and  y  *  -jj , 


there  results 


[36] 


For  two  different  pressure  gradients  at  the  same  Reynolds  number  R„  and  conse 


quently  at  the  same  value  of  C  and  n  , 
the  following  ratio  can  be  formed 


If  the  reference  condition  1b  that 
for  zero  pressure  gradient,  this 
equation  becomes 


where  the  subscript  zero  refers  to 
zero  pressure  gradient.  For  a  power- 
law  velocity  profile  of  a  flat  plate 
with  zero  pressure  gradient.  Equation 
[14]  becomes 


Figure  7  -  Regions  of  Flow  Within  the 
Turbulent  Boundary  Layer  in  a 
Pressure  Gradient 
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Figure  8  -  Logarithmic  Plot  of  Velocity  Profiles  in  Pressure  Gradients  a 
Illustrating  "The  Law  of  the  Wall"  (Ludwleg  and  Tlllmann,  Reference  18) 


Then  [38]  reduces  to 


[40] 


Comparing  the  exponents  of  y  in  [34]  and  [40],  it  Is  seen  that  the  exponent 
obtained  by  Ludwleg  and  Tlllmann  has  a  constant  value  of  1 .703  while  that  of 
the  equation  Just  derived  depends  on  H0,  which  varies  with  Reynolds  number. 
Prom  Figure  11  showing  the  variation  of  Hp  with  R# ,  the  exponent  4/H#+l  is 
found  equal  to  1.639  for  R(  «  1 ,5  x  10s,  and  1 .800  for  R#  «  4  x  104.  These 
values  compere  favorably  with  the  value  1 .705  of  Ludwleg  and  Tlllmann. 


1.4 


Since  from  Equation  [9] 

y  «  Vi  -  n  [41] 

y  may  be  expressed  as  a  function  of  H  using  [15] • 


r 


[42] 


and 


o  Lh0(h.+i)J 


[43] 


(y  \§J+\ 

j-J  0  with  H  and  H0  is  shown  in  Pigure  9.  Since  this  fac¬ 
tor  which  represents  the  effect  of  pressure  gradient  on  the  skin  friction  di¬ 
minishes  asymptotically  towards  zero,  with  increasing  H,  it  can  only  describe 
flow  conditions  approaching  separation  and  not  conditions  right  at  the  separa¬ 
tion  point  where  the  wall  shearing  stress  is  zero. 
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LOCAL- SKIN-FRICTION  COEFFICIENTS  FOR  FLAT  PLATES 
AT  ZERO  PRESSURE  GRADIENT 

Since. the  evaluation  of  the  local  skin  friction  in  a  pressure  gradi¬ 
ent  depends  on  a  knowledge  of  its  value,  in  zero  pressure  gradient,  as  seen  in 
Equation  [40],  a  review  of  available  information  on  the  drag  of  flat  plates 
should  now  be  considered.  A  smooth  flat  plate  moving  parallel  to  the  direc¬ 
tion  of  motion  represents  an  excellent  example  of  boundary-layer  flow  in  a 
zero  pressure  gradient.  The  analytic  study  of  the  drag  of  flat  plates  has 
proceeded  on  a  semi-empirical  basis  by  such  investigators  as  Prandtl  and 
von  Karman.1’2  From  various  theoretical  and  empirical  considerations  a  func¬ 
tional  form  is  assumed  for  the  velocity  profile  involving  both  the  velocity 
and  the  local  skin  friction.  Then  an  integration  over  the  length  of  the  plate 
leads  to  a  general  expression  for  the  drag,  containing  coefficients  to  be 
numerically  evaluated  from  test  data.  Landweber21  has  made  a  critical  review 
of  such  methods. 

For  the  boundary-layer  relations  considered  in  this  report  it  is 
necessary  to  express  the  frictional  resistance  of  flat  plates  in  terms  of  a 
local-skln-frictlon  coefficient  as  a  function  of  a  local  Reynolds  number  or 


^  -  f(R,)  [44] 

P  U2 


where  R#  *  Various  Investigators  have  developed  expressions  of  this  type 
from  drag  coefficient  formulas  appearing  in  the  literature. 

Tetervln22  obtained  the  following  local-skin-frlction  formula  for 
flat  plates  from  the  universal  resistance  law  for  pipes 


Since  rw  //>U2  appears  also  on  the  right  hand  side,  it  is  an  unwieldly  equation 
to  use  numerically.  Furthermore,  as  it  is  based  on  pipe-resistance  data,  it 
applies  only  approximately  to  flat  plates. 

Squire  and  Young11,  by  a  combination  of  the  von  Ka'rman  asymptotic 
formula  and  the  Prandtl  and  Schlichtlng  drag  law,  derived  the  following  ex¬ 
pression  for  the  local-skin-frictlon  coefficient 

~$m  [5-890  logj (4.075  R#>] 


[23] 
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The  Prandtl  and  Schlichting  drag  law  has  been  criticized  for  its  use  of  pipe 
velocity  profiles  in  its  derivation.3 

Palkner,13  by  a  study  of  various  test  data  for  flat  plates,  proposed 
the  following  simple  power -law  formula 

li*.  ■  0-^5|S  [25] 

oV2  R,1'* 

Like  all  power-law  formulas,  this  equation  is  valid  only  for  a  limited  range 
of  Reynolds  numbers. 

Ludwieg  and  Tillmann18  derived  the  following  equation  from  the  test 
data  of  Schultz-Grunow23 


0.0167 

Uog10R,)1M8 


[46] 


for  2  x  103  <R#  <7  x  105. 

A  local-skln-friction  formula  will  now  be  derived  from  the  well- 
known  Schoenherr  (K^rman-Schoenherr)  frictional-resistance  formula  for  flat 
plates.  Prom  an  extensive  collection  of  test  data  for  the  frictional  resist¬ 
ance  of  flat  plates  covering  a  wide  range  of  Reynolds  numbers,  Schoenherr24 
supplied  numerical  values  to  the  coefficients  of  the  von  Ktfrman  asymptotic 
drag-law  formula.  The  Schoenherr  formula  has  become  the  baslB  of  frictional- 
resistance  calculations  for  the  prediction  of  full-scale  horsepower  from  tests 
on  models  of  ships.  The  Schoenherr  formula  is 


-^T  -  4.13  log jj  (RxCf ) 


[47] 


where  Rx  •  Here 


[48] 


where  D  is  the  drag  and  S  is  the  total  surface  area.  Hence  for  a  flat  plate 
of  width  W  and  length  x  the  drag  of  one  side  is 


D  -  Cf  ^>U2x  W 


[49] 
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Foi  a  flat  plate  with  zero  pressure  gradient  the  von  Krfrman  momentum  equation 
reduces  to 


41.  T*o 

W 


.[50] 


or 


or 


where  by  definition 


Hence 


9 


r  dx 

"o 


9  - 


1  D 


x 

b  -  w  /  r  ax 
0  0 


D  -  pU*«W 


(51) 

[52] 


[53] 


Finally  from  Equations  [49]  and  [53] 


C 


f 


[54] 


The  substitution  of  the  relations  for  Cf  from  [54]  and  for  R^  into  the 
Schoenherr  formula,  [47],  produces 


[55] 


Differentiation  with  respect  to  x  gives 
d« 


(9\l  ^ 

r  i 

8.26  VT 

lo8io(^r)*  +  loBu* 

[56] 


Substituting  for  from  [55]  and  using  [50]  finally  results  in 


Tw0  m  _ 0.01466 _ 

pU*  logw(2R#)[^  log10(2R#)+0.4343] 


[57] 


With  the  exception  of  Falkner's  equation,  all  the  local-skin- 
friction  formulas  discussed  in  this  section  are  plotted  together  in  Figure  10 
for  comparison.  It  is  seen  that  there  are  no  appreciable  numerical  differ¬ 
ences  among  than.  For  convenient  reference  the  flat-plate  formulas  are  listed 
in  Table  2 . 


TABLE  2 


Summary  of  Formulas  for  Local- Skin-Friction  Coefficients 

of  Flat  Plates 


Investigator 

Date 

Origin 

Formula  for  Local-Skln-Frlctlon  Coefficient 

Remarks 

Tetarvln 
Reference  22 

19A4 

Universal  Pipe  Re¬ 
sistance  Lav 

m 

Based  on  pipe  re¬ 
sistance  data 

mm 

Squire  and 
Young 

Reference  11 

1937 

von  Kirman  Asymp¬ 
totic  Formula  and 
Prandtl  and  , 
Schlichting  Drag 

Law 

p[}>  *  (5.890  iogu4.075  H,)  [23] 

Fallen  er 
Reference  13 

19**3 

Experimental  Data 

VjLOS&ii  1251 

FU*  R, 

Limited  renge  of 
Reynolds  number 

Ludwleg  and 
Tlllaann 
Reference  18 

1949 

Schults-Orunow 

Data 

\  .  -,<Lg]6L-  [46] 

«*  (logMR,) 

OranvlUe 
(This  report) 

1951 

Schoenherr  Formula 
(von  Xdrmin  Asymp¬ 
totic  Formula) 

- ;  (57) 

W  1H„  UK,  Jfylogj,  (2R,  )+0 .4343] 

In  addition  to  the  local-skin-frlction  coefficient,  the  variation  of 
the  shape  parameter  of  the  velocity  profiles  Hfl  of  flat  plates  with  Reynolds 
number  at  zero  pressure  gradient,  is  required  in  Equation  [40]  to  evaluate  the 
local-skin-friction  coefficient  in  a  pressure  gradient.  From  the  measurements 
of  Schultz-Grunow,23  and  from  unpublished  British  data,  Tetervln  and  Lin15 
formulated  the  following  empirical  expression 

logX0H0  -  0.5990  -  0.1980  logl#R,  +  0.0189  (logwR,)2  [58] 

for  1 .5  x  10*  <  R#  <  1  x  105.  Figure  11,  where  Equation  [58]  is  plotted, 
shows  the  decrease  in  Hg  with  increasing  Reynolds  number  R#. 
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Figure  10  -  Local-Skin-Frictlon  Coefficient  t, /pU*  as  Function  of 
Local  Reynolds  Number  Rt  for  a  Flat  Plate 
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INTEGRAL  OP  THE  SHEARING -STRESS  DISTRIBUTION 
ACROSS  THE  BOUNDARY  LAYER 


In  order  for  the  moment-of -momentum  equation  [27]  to  be  used  numer¬ 
ically  a  relationship  is  required  which  describes  the  effect  of  pressure  gradi¬ 
ents  on  the  integral  of  the  shearing-stress  distribution  across  the  boundary 
layer j  Although  almost  no  study,  either  theoretical  or  empirical, 

has  been  made  of  the  Integral  of  the  shearing-stress  distribution,  some  stud¬ 
ies  have  been  made  of  the  shearing-stress  distribution  Itself. 

The  hot-wire  measurements  of  Schubauer  and  Klebanoff,85  shown  in 
Figure  12,  demonstrate  the  large  changes  produced  in  the  shape  of  the  shearing- 
stress  profile  by  an  adverse  pressure  gradient.  The  shearing-stress  curves 
have  a  positive  slope  at  the  wall  determined  by  the  positive  pressure  gradient; 
they  then  come  to  a  peak,  and  finally  drop  to  zero  at  the  outer  edge  of  the 
boundary  layer.  Several  attempts  have  been  made  to  derive  analytical  expres¬ 
sions  for  the  general  shearing-stress  distribution;  these  are  described  below. 

Fediaevsky,2*  following  the  Pohlhausen  method  for  laminar  flow, 
tried  to  fit  a  polynomial  expression  for  the  shearing-Btress  distribution  that 
would  satisfy  the  boundary  conditions  at  the  inner  and  outer  edges  of  the 
boundary  layer.  The  coefficients  Aj  of  the  following  polynomial  are  to  be 
evaluated: 


< 

f  ■  z  *  .(f/ 

rw  J-0 

[59] 

The  boundary 

conditions  imposed  by  Fedi&evsky  are: 

a. 

At 

y  « 

o, 

JL-  1 

Tw 

by  definition, 

b. 

At 

y  * 

o, 

6rB 
dy  3x 

from  the  boundary-layer 
equation,  Equation  [1], 

,c. 

At 

y  - 

a. 

T  s  0 

from  the  definition  of 
boundary-layer  thickness , 

d. 

At 

y  ■ 

a. 

lr«  o 

By 

from  the  assumption  that 
the  derivative  of  the  total 
head  of  the  outer  streamline 
is  continuous, 

e. 

At 

y  - 

o, 

ft*  0 

ay1 

from  the  derivative  of  the 
boundary-layer  equation , 
Equation  [1  j  and  a  p/ay  ■  0. 
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Substituting  only  the  first  three  boundary  conditions  a,  b,  and  c  into  Equa¬ 
tion  [59]  and  its  derivatives  to  form  a  system  of  simultaneous  equations  for 
evaluating  the  coefficients,  results  in  the  following  quadratic  equation  for 
the  shearing-stress  distribution 

—  ■  1  +  «(^r)  -  0  +  [60] 

rw 

where  a  •  ~  ■'Jf,  a  pressure-gradient  parameter.  Using  all  five  boundary  con¬ 
ditions  yields  the  following  quadratic  equation  for  the  shearing-stress 
distribution 


rw  “  1 


+  a 


(-f)  -  <*  ♦  3.)  (-f  )*  +  (3  ♦  2.)  (-f  )*  1611 


Schubauer  and  Klebanoff*®  found  indifferent  correlation  between  their  experi¬ 
mental  results  and  either  of  Pediaevsky's  equations.  Agreement  is  good  at  the 

beginning  of  the  adverse  pressure  gradi¬ 
ent  and  near  the  point  of  separation 
biit  is  extremely  poor  at  the  inter¬ 
mediate  stations. 

In  an  attempt  to  improve  on 
the  Fedlaevsky  formulations,  Ross  and 
Robertson® 7  considered  the  shearing- 
stress  distribution  to  depend  on  its 
previous  history  of  development  along 
the  boundary  layer.  Pediaevsky's  bound¬ 
ary  condition  (d)  is  altered  to 
at  y  ■  4,  in  order  to  obtain  a  constant 
slope  at  the  outer  edge  of  the  boundary 
layer  travelling  downstream  from  1,  the 
initial  position  of  the  adverse  pres¬ 
sure  gradient.  This  requirement  makes 
the  shearing-stress  distribution  depend 
on  conditions  upstream.  Using  a  modi¬ 
fied  polynomial  expression  for  r  and 
incorporating  the  revised  boundary  con¬ 
dition  (d),  together  with  the  other 
four  Fedlaevsky  conditions,  Ross  and 
Robertson  derived  the  following 
expression: 


1 

a 

91 

9| 

a 

m 

m 

n 

m 

a 

m 

m 

B 

a 

y  in 

Figure  12  -  Shearing- Stress 
Distributions  in  an  Adverse 
Pressure  Oradlent  as 
Measured  by  Schubsuer  and 
Klebanoff,  Reference  25, 
At  Various  Stations 


where 


r 


[62] 


V 


\ 


( 
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4’*#)-  ,•♦/*)(!  -f)f  ♦<-♦-»('  -f)' 


?  +  l 


and 


9  .  ?«  ±  l  ±  1 
*  -  1 


A  comparison  of  the  above  expression  with  the  experimental  findings  of 
Schubauer  and  Xlebanoff  shows  a  somewhat  better  agreement  than  the  Fedlaevsky 
relations . 

The  integral  of  the  shearing-stress  distribution  will  now  be  studied 
by  integrating  the  analytic  relations  just  described  for  the  shearing-stress 
distribution  and  by  analyzing  the  test  data  of  Schubauer  and  Klebanoff .  Let 


1  *  (n  )  (63 I 

By  elementary  integration  the  following  expressions  are  derived:  Por  the 
Fedlaevsky  3-condltlon  polynomial,  Equation  [60], 

I  -  0.67  +  0.17a  [64] 

for  the  Fedlaevaky  3-condltlon  polynomial ,  Equation  [6l ] , 

1  -  0.60  +  0.13a  [63] 

and  for  the  Ross  and  Robertson  equation,  Equation  [62] 

,  1  (66] 

c  '  2a  +  -  1  • 


Significant  results  are  obtained  from  the  data  of  Schubauer  and 
Klebanoff  when  the  Integral  of  the  shearing-stress  distribution  Is  non- 
dlmenslonallzed  as  follows 


[68] 
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Values  of  are  computed  directly  from  the  measurements  of  Schubauer  and 
Klebanoff  by  substituting  data  into  Equation  [67]  and  performing  the  indicated 
Integrations  numerically  by  Simpson's  rule.  Por  comparison,  values  of  #  are 
obtained  from  the  analytical  formulations  of  I  based  on  the  aethods  of 
Pediaevsky  and  of  Ross  and  Robertson  by  inserting  values  of  a  and  fi  from  the 
experimental  data  into  Equations  [64],  [65],  and  [66].  As  seen  in  Figure  13, 
in  which  tf>  is  plotted  for  various  stations  along  the  boundary  layer,  both  of 
the  Pediaevsky  formulations  show  poor  correlation  with  experimental  points 
while  the  Ross  and  Robertson  formulation  shows  only  fair  agreement.  It  is 
observed  that  except  for  experimental  error  the  value  of  $  is  practically  con* 
stant  for  different  stations  along  the  boundary  layer.  This  suggests  that 
using  a  value  of  ft  independent  of  the  effects  of  pressure  gradient  ought  to  be 
sufficiently  accurate  for  most  technical  problems  encountered. 


Before  deciding  on  values  of  f  for  pressure  gradients  in  general, 
it  will  prove  fruitful  to  investigate  the  values  of  ^  in  a  zero  pressure  gradi¬ 
ent.  Proa  Equation  [68]  i>  for  a  zero  pressure  gradient  or  #0  may  be  written 


as 
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Pros  [12],  where  a  power  law  for  the  velocity  profile  Is  assumed, 


Then 


t  H.  +  1 


[70] 


■  §)(£)’ 


[71] 


The  value  of  I#  Is  now  derived  analytically  from  the  momentum  and 
from  the  moment -of -momentum  equations  for  flat  plates  with  zero  pressure  gradl 
ent  ( ^  -  0)  In  the  following  manner.  The  moment -of -momentum  equation,  [27], 
reduces  to 


d 


<Ho 


i)[h0 


(H„ 


and  the  momentum  equation  to 


d  9 
dx 


[72] 


[50] 


Combining  both  equations  gives 


9Wm 

•  '>[“. 

+  ’)  1.] 

da 

dx 

or 

d«o 

(H*  -  D[i 

i,  -  <«,  + ')  i, 

1 

da  * 

9 

In  terms  of  R#  -  ^r,  [74] 

becomes 

dH, 

(H?  -  1) 

H0  -  (H0  +  1)  ] 

dR# 

R, 

[73] 

[74] 


[75] 


Differentiation  of  the  empirical  equation  of  H0,  Equation  [58],  with  respect 
to  R(  produces 

dH0  (-0.1980  +  0.0578  log10R,)  Hp 

dS7  T, 


[76] 
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Finally  I#  is  evaluated  by  equating  [75]  and  [76].  Then 


V1 


1  + 


0.1980  -  0.0378  log1#R# 


Ho  -  1 


[77] 


R,  may  be  eliminated  by  use  of  Equation  [58]  to  give 

H. 


H0+1 


1  + 


0.0378  y$2.9  log10H0  -  47T8 


178] 


The  small  variation  of  I#  with  Reynolds  number  R#  is  shown  in  Figure 
14.  For  instance,  at  R#  =  T500,  I0  =  0.633  and  at  R#  =  100,000,  I  ..  0.554. 
These  values  appear  reasonable  when  compared  with  those  based  on  the  Fediaevsky 
formulations:  From  Equation  [64]  for  three  conditions  I  =  0.67  (a  =  0  for 
zero  pressure  gradient)  and  frojn  Equation  [65]  for  five  conditions  I  «  0.60. 


TJ 


4? 


o 


V 


i 


Finally  from  Equations  [71]  and  [78],  is  written 

0 . 0378 /52.9  log10 H0  -  4.18 


Using  the  value  of  /p\jz  from  [57],  derived  from  the  Schoenherr 
formula,  and  the  value  of  Hp  from  [58],  the  variation  of  with  Reynolds 
number  is  shqwn  graphically  in  Figure  14.  It  is  seen  that  diminishes  grad¬ 
ually  with  Reynolds  number. 
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When  the  values  of  0#,  based  on  the  same  Reynolds  number  as  those 
appearing  In  the  experiments  of  Schubauer  and  Klebanoff ,  are  plotted  In  Figure 
13  for  comparison,  they  produce  a  curve  of  almost  constant  value  which  Is  In 
excellent  agreement  with  the  experimental  points  of  Based  on  this  evidence 
it  will  be  reasonable  to  conclude  that  at  the  same  Reynolds  number  Is  a 
close  approximation  to  0  for  other  boundary  layers  In  pressure  gradients  sim¬ 
ilar  to  the  boundary  layer  of  Schubauer  and  Klebanoff.  Further  substantiation 
depends  on  additional  experimentation,  especially  at  high  Reynolds  numbers. 

The  integral  of  the  shearing-stress  distribution  as  required  by  the 
moment -of -momentum  equation,  Equation  [27],  is  expressed  in  terms  of  ^  as 
follows : 


For  power-law  velocity  profiles 


H  -  1 

V"  FTT 


[12] 


and  with  the  assumption  -  *#,  where  Is  defined  In  Equation  [71  ] ,  the 
Integral  of  the  shearing- stress  distribution  finally  becomes 


This  completes  the  shearing-stress  relations  required  for  transform 
lng  the  moment -of -momentum  equation  Into  an  auxiliary  equation  for  character¬ 
izing  the  shape  of  the  velocity  profiles. 


MODIFIED  MOMENTUM  AND  MOMENT-OF-MOMENTUM  EQUATIONS 


The  substitution  of  the  relation  for  local  skin  friction  In  a  pres¬ 
sure  gradient,  Equation  [40],  Into  the  von  Karmdn  momentum  equations,  [6a]  and 
[6b],  modifies  them  as  follows:  For  two-dimensional  flows, 

(«.) 

and  for  axlsymmetric  flows  (where  *4 crw), 


dx 


-  <»  ♦  *>  4  3H£) 


4 

FS+T 


[82b] 
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The  Insertion  of  the  relations  for  local  skin  friction  and  for  the  integral 
of  the  shearing-stress  distribution,  [40]  and  [8l],  into  the  moment-of- 
momentum  equation  [27],  converts  it  into  the  following  equation,  which  is  ap 
plicable  to  the  calculation  of  both  two-dimensional  and  axlsymmetrie  flows 
where  ^«rw) 

(*r- 

where  y,  and  Io  are  defined  in  Equations  [42],  [43],  and  [78]  respectively, 
and  the  subscript  0  refers  to  flat-plate  values  at  the  same  Reynolds  number  R#. 
This  equation  is  the  desired  auxiliary  equation  to  be  solved  simultaneously 
with  one  of  the  modified  von  Karman  momentum  equations,  [82a]  or  [82b],  for 
calculating  the  growth  of  boundary  layers  in  pressure  gradients. 

The  modified  moment -of -momentum  equation  of  this  report  represents 
an  auxiliary  equation  whose  form  has  a  theoretical  basis  and  whose  coeffi¬ 
cients  have  been  evaluated  indirectly  from  flat-plate  data.  On  the  other  hand 
the  various  auxiliary  equations  described  in  this  report  are  almost  wholly 
empirical  in  origin.  A  comparison  of  the  auxiliary  equations  may  be  made  by 
writing  them  in  the  following  generalized  form 

»§-■-*  (vf)-B  w 


«%-  (n..n 


The  expression  for  the  variation  of  rj  in  the  Oruschwitz  or  Kehl  method  is  con¬ 
verted  into  the  above  form  by  using  the  power-law  formula,  [15],  for  the  rela¬ 
tion  between  n  and  H.  The  expressions  for  the  coefficients,  A  and  B,  are 
listed  in  Table  3  for  the  various  auxiliary  equations.  It  is  seen  that  the 
coefficients  for  the  Kehl,  Garner  and  the  moment -of -momentum  equations  are 
functionally  similar  with  A  -  f(H)  and  B  -  f(H,  R# ) .  The  von  Doenhoff  and 
Tetervin  equation,  nowever,  has  the  Reynolds-number  effect  in  coefficient  A 
Instead  of  coefficient  B.  As  expected,  none  of  the  coefficients  for  the 
Gruschwitz  method  exhibits  a  Reynolds-number  effect. 

The  coefficients  A  and  B  are  plotted  against  H  in  Pigures  15,  16, 
and  17  for  two  different  values  of  R#.  It  is  observed  that  although  all  the 
curves  have  the  same  general  shape  and  trend,  there  are  marked  differences  in 
the  magnitude  of  the  ordinates.  The  exponential  form  of  the  coefficients  of 
the  Garner  and  of  the  von  Doenhoff  and  Tetervin  equations  is  responsible  for 
the  very  large  increase  in  value  of  the  coefficients  at  the  higher  values  of 
H.  The  curves  for  the  coefficients  of  the  modified  moment -of -momentum  equa¬ 
tion  of  this  report  maintain  average  positions  among  the  other  curves. 
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Figure  17  -  Coefficient  B  of  General  Auxiliary  Equation  Compared  for 
Various  Formulations  at  Reynolds  Number  R(  «  1  x  105 

Although  this  fact  does  not  conclusively  prove  the  validity  of  the  modified 
moment -of -momentum  equation,  It  does  indicate  that— as  such— the  modified 
moment -of -momentum  equation  has  a  proper  range  of  values  for  an  auxiliary 
equation  giving  the  variation  of  H.  There  is  a  need  for  a  larger  body  of  re¬ 
liable  measurements  of  turbulent  boundary  layers  under  a  variety  of  conditions 
and  for  increased  emphasis  on  measurements  of  the  shearing  stresses  both  at 
the  surface  and  throughout  the  boundary  layer. 


SUMMARY  OF  METHOD  OF  CALCULATING  TURBULENT  BOUNDARY  LAYERS 
IN  A  PRESSURE  GRADIENT 


For  calculating  turbulent  boundary  layers  in  a  pressure  gradient  the 
following  relations  are  involved:.  The  momentum  equation  for'  two-dimensional 
flows, 


-|h  +  j)T35 


de 

dx 


[82a] 


or  the  momentum  equation  for  axisymmetric  flow  (i«rw), 


d(ri!) .  .(a  ♦  2) :«! i"+(|)vr  Is. 

U  dx  'V  " 


dx 


[82b] 


pU‘ 


and  the  moment -of -momentum  equation  which  holds  for  both  two-dimensional  and 
axisymmetric  flows  (<<<r  ), 


•*  -  <«  -  ,)|^}  i, 

-  o  0 

where 


(Hft+1 ) 


^  [83] 
“J  pU2 


r  ■  [sfirn]^ 

[42] 

and 

f  H0-1  1  «“ 

y«“  k(H0+1)J 

[43] 

The  variation  of  (Vyo)*o+1  with  H  and  H0  ls  sh°wn  ln  Figure  9.  The  varia¬ 
tion  of  H0  with  R4,  shown  ln  Figure  11,  is  given  by  Reference  15  as 


log10H0  =  0.5990  -  0.1980  log10R,  +  O.Oi89(log10R,)s 


The  relation  between  IQ  and  H0'  is 


I  =  — — 

'  o  H„+1 


0  L 


1  + 


0.0378  F52.9  iog10H0-4.i8 


-  1 


[58] 


[78] 


Figure  14  shows  I„  as  a  function  of  R, .  The  variation  of  rw  /pU2  with  R#, 
Figure  10,  for  the  Schoenherr  formula  is 


0.01466 


log10(2R,)[y  1og10(2R,)+0.4343] 


[57] 


The  variation  of  U  and  dU/dx  with  x  will  be  specified  for  a  particular  problem 
and  may  be  obtained  either  by  pressure  measurements  or  potential-flow 
calculations . 

The  momentum  and  moment -of -momentum  equations  are  solved  for  0(x) 
and  H(x)  as  a  pair  of  simultaneous  differential  equations,  usually  by  numer¬ 
ical  methods  Involving  step-wise  Integration.  The  Initial  conditions  for  B(x) 
and  H(x)  are  Imposed  by  the  physical  conditions  of  the  problem  such  as  the 
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transition  point  from  laminar  to  turbulent  flow.  It  Is  customary  to  use  the 
value  of  0  of  the  laminar  boundary  layer  at  the  transition  point  as  the  Ini¬ 
tial  value  of  0  for  the  turbulent  boundary  layer.  Owing  to  the  lack  of  reli¬ 
able  data  on  transition,  the  Initial  value  of  H  may  be  taken  as  that  for  a 
flat  plate  at  the  -same  Reynolds  number  R(. 

The  evaluation  of  $  and  H  is  required  in  many  technical  hydro- 
dynamic  and  aerodynamic  applications  involving  calculations  of  drag  and  of 
velocity  profiles  of  boundary  layers.  According  to  the  following  formula  of 
Squire  and  Young11  for  the  profile  drag  of  airfoils,  0  and  H  must  be  deter¬ 
mined  at  the  trailing  edge  H  +5 

/ 

»  -  pvi'STJ  [85] 

where  the  subscript  t  refers  to  conditions  at  the  trailing  edge.  If  a  power- 
law  distribution  is  assumed,  the  values  of  0  and  H  at  any  station  along  the 
boundary  layer  specify  the  shape  and  thickness  of  the  velocity  profile,  as 
shown  by  Equation  [l6j, 


H-  1 


A  knowledge  of  the  characteristics  of  the  velocity  profiles  of  boundary  layers 
is  required  for  the  design  of  control  surfaces  and  propellers  on  forms  moving 
in  fluids  and  of  the  various  conduit  sections  of  water  and  wind  tunnels. 

Large  values  of  H  indicate  the  danger  of  flow  separation  for  cases  involving 
adverse  pressure  gradients. 
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